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Abstract The chromospheric response to the input of flare energy is marked by ex¬ 
tended extreme ultraviolet (EUV) ribbons and hard X-ray (HXR) footpoints. These are 
usually explained as the result of heating and bremsstrahlung emission from acceler¬ 
ated electrons colliding in the dense chromospheric plasma. We present evidence of 
impulsive heating of flare ribbons above 10 MK in a two-ribbon flare. We analyse the 
impulsive phase of SOL2013-11-09T06:38, a C2.6 class event using data from Atmo¬ 
spheric Imaging Assembly (AIA) on board of Solar Dynamics Observatory (SDO) and 
the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) to derive the 
temperature, emission measure and differential emission measure of the flaring regions 
and investigate the evolution of the plasma in the flaring ribbons. The ribbons were 
visible at all SDO/AIA EUV/UV wavelengths, in particular, at 94 and 131 A filters, 
sensitive to temperatures of 8 MK and 12 MK. Time evolution of the emission measure 
of the plasma above 10 MK at the ribbons has a peak near the HXR peak time. The 
presence of hot plasma in the lower atmosphere is further confirmed by RHESSI imag¬ 
ing spectroscopy analysis, which shows resolved sources at 11-13 MK associated with 
at least one ribbon. We found that collisional beam heating can only marginally explain 
the necessary power to heat the 10 MK plasma at the ribbons. 
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1. Introduction 


The impulsive phase of solar flares is characterised by intense emission of microwaves 
and hard X-rays (HXR), showing the presence of non-thermal electrons in the corona 
and chromosphere, respectively. These electrons are believed to be accelerated by the 
magnetic energy release mechanisms in the coron a (e.g. Z hark ova et al. llOl ll) although 
this interpretation has bee n contested recently dFletcher and Hudsonll2008l : lBrown et al. 


2009l : IVaradv etal. 112014b . Strong ultraviolet (UV) and extreme ultraviolet (EUV) emis- 


sion is also commonly observed in association with the non-thermal emission (e.g. 


Hinteregger and HalJll969l : lEmslie and Novesl[l978l : iHoran. Kreplin. and Fritz!Il982l : 
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Alexander and Covnei]l2006l : ICovner and Alexanded 12009). indicating fast heating at 


the transition region and chromosphere. This heating is usually explained by the en¬ 
ergy depo sition of the accelerated ele c trons hitting the hig h-density plasma at these 
locations dBrown|[l973l : lHudsonlll972l : iBrown et <^Z.I Il990h . The increase of pressure 
due to heating drives the plasma upwards into the coronal loops, a process termed 
chromospheric evaporation. The coronal loops filled with hot plasma will be bright 
in soft X-rays_(SXR) and EUV, peaking at temperatures between 8 and 40 MK ( e.g. 
Rvan et <2/1 [2012b . The maximum temperature is usually reached after the impulsive 


phase, i.e. after most of the released energy is deposited in the ambient plasma. The hot 
coronal loops eventually cool by conduction, sending energy to the lower atmosphere, 
followed by_a long radiative cooling phase (e.g. ISves tkalll 987th . 

McTieman et al. dl993h . using Yohkoh Soft X-ray Telescope (SXT) images identi- 
fied impulsive SXR emission associated with a HXR footpoint. Its duration was less 
than one minute, reaching its maximum about 20 seconds before an associated HXR 
sub-burst was detected by the low energy channel, and 40 seconds before the main 
HXR burst at the same location. This SXR burst corresponded to a sharp temperature¬ 
time profile, risin g fro m ^_7 MK to about « 10 MK in « 30 seconds, followed by 
a fast cooling. [Hudson et al\ ( 1994 ) reported similar observations for the limb event 
SOL1992-01-26, where the footpoint SX R peak time matched the HX R peak time 
within the instrument’s temporal resolution. llMrozek and Tomczak ( 2004 ) studied SXR 
impulsive brightenings in footpoints in 46 Yohkok events, and concluded that these 
can be explained as result of collisional heating mainly by non-thermal electrons with 
relatively low energies_ 

More recently, [Graham et all d2013l) presented the emission measure distributions 
(EMD) as a function of temperature for the footpoints of six flares, derived from EUV 
spectroscopic observations, and showed that the footpoint EMDs had a peak around 8 
MK, indicating the presence of high temperature plasma within the footpoints. Evidence 
of hot (T > 1 MK) and dense (n e > 10 10 cm -3 ) plasma at fl are fo ot points d uring the 
impulsive phase was also presented by several other authors dMilligan and Dennisl2009l : 
Watanabe et al. Eoiol: Del Zanna et al. 201 ll: iMilliganl 201 ll: Graham. Fletcher, and 


Hannah 201 ll : iFletcher et al\ 2013b . A comprehens ive review of E UV spectroscopic 
observations of hot footpoints has been presented by Milligan) ( 2015 ) 


While the Hino del Extreme ultraviolet Imaging Spectrometer (EIS: Culhane et al 


20071) provides great spectroscopic information on the plasma across a wide temper¬ 


ature range, the slit-raster observational technique does not allow evaluation of the 
plasma evolution at the same spatial location with high t empor al cadence. Imaging 
from the Atmospheric Imager Assembly (AIA: iLemen et al. 2012 ). on board the Solar 
Dynamics Observatory (SDO), and methods to recover the EMDs allow us to track the 
plasma evolution at different regions of flares simultaneously. We present the analysis 
of the SOL2013-11-09 flare ribbons, recovering the temperature and emission measure 
and show evidence for a sudden heating to « 10 MK temperatures of the ribbon plasma 
during the impulsive phase. 


2. Observational Data of SOL2013-11-09 

The SOL2013-11-09 flare occurred in active region NOAA 11890 near disc centre, and 
was characterised by two flaring ribbons plus a bright compact source between them 
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during the impulsive phase. For this study we used data from SDO/AIA. AIA has nine 
filters, or passbands, in theEUV andUV range: 94, 131, 171, 193, 211, 304, 335, 1600 
and 1700 A filters, taking full Sun images with a pixel resolution of 0.6 arcsecs and 12 
seconds cadence for the EUV passbands and 24 seconds for the UV passbands. We also 
used data from the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI: 


Lin et a/.ll2002l) . The RHESSI imager consists of nine bi-grid rotating modulation col¬ 


limators (RMCs), with a rotation period of about 4 seconds. The data from individual 
RMC are combined by the imaging algorithms to cover different spatial scales, from 
2.3 arcsecs (RMC 1) up to 180 arcsecs (RMC 9). To avoid saturation of the detectors at 
times with high rate of incident photons, two-stage (“thin” and “thick”) attenuators can 
be automatically employed. No attenuators (state “A0”) were in during this event, a rare 
occurrence that gives us the opportunity to investigate the evolution of the low-energy 
X-rays during the impulsive phase without gaps and with full count rates. Complemen¬ 
tary SXR data from the Soft X-ray Sensor on board of the Geostationary Operational 
Environmental Satellite (GOES) were also employed. GOES observes the Sun as a star 
in two broadband SXR channels, 1-8 A and 0.5-4 A, which can be used to estimate the 
temperature, T, and emission measure, EM , o f the flaring plasma ( Thomas. Crannell, 
and Starr ll985l [White, Thomas, and Schwartzll20051) . "" 

In Figure |T] we summarise some aspects of the event. Figure [T^ shows RHESSI 
HXR count rates at several energy bands, where the high-energy bands (9-12, 12-25 
and 25-50 keV) show a clear impulsive peak with a maximum at 06:25:46 UT, and are 
associated with bremsstrahlung emission from non-thermal electrons, as we discuss in 
Section® A second, weaker and less impulsive peak is seen around 06:27:20 UT. The 
lower energy bands (below 9 keV), associated with thermal bremsstrahlung, show a 
gradual rise starting slightly after 06:22 UT, and also a pronounced peak simultaneous 
with the high-energy channels. Both GOES channels, 1-8 A and 0.5-4 A (Figure 
Ht) have a steep rise, with a “shoulder” at the same time as the HXR peak. Spatially 
integrated EUV emission was obtained by summing the pixels of AIA images taking the 
whole flare region and are shown in Figure®:. The 171 A emission is well associated 
with the high-energy HXR, with a main impulsive peak around 06:25:46 UT and a 
secondary peak around 06:27:20 UT. Similar characteristics are observed at 193, 211, 
304, 1600, and 1700 A. Emission in the 94 and 131 A passbands is similar to the GOES 
SXR channels: a sharp rise matching the HXR peak, followed by a more gradual rise. 
These AIA passbands are sensitive to emission from hot plasma: the main contribution 
to the 94 A passband is the Fe xviii 93.93 A emission line, formed at logT = 6.8 
while the dominant contributio n to the AIA 131 A passband during flare s comes from 
Fe xxi 128.75 A (log T = 7.1) dO’Dwver et a/.ll2010l : lBoerner et fl/.M2014l) . These main 


contributions were verified using from the EUV Variability Experiment (EVE: Woods 


et < 


al. 120121) data for this event by Simoes, Graham, and Fletcherl d2015l) . 


3. SDO/AIA Data Analysis 

3.1. Impulsive EUV Emission from Flaring Ribbons 

We will now focus our analysis on the impulsive phase, 06:22 to 06:30 UT, when the 
impulsive ribbon emission is observed. AIA images (see Figure [2] and [3]) reveal two 
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Figure 1 . SOL2013-11-09 flare overview, (a) RHESSI HXR count rates (energy bands indicated in the 
figure), (b) GOES flux at 1-8 and 0.5-4 A channels, (c) Spatially integrated emission of AIA passbands 
94, 131, 171 A. Bottom: (d) Temperature, T, and emission measure, EM, determined using from GOES and 
RHESSI. In all frames, the vertical dotted line indicates the HXR peak. 


ribbons as the brightest features in all AIA EUV/UV passbands during the impulsive 
phase of the event. The ribbons are located in regions with opposite magnetic field 
polarity, as can be seen in Figure ® which shows the line-of-sight (LOS) photospheric 
mag netic field obtained bv SDO IHelioseismic and Magnetic Imazer (HMI: Scherrer 


et al. 20121) . Coronal loops become the dominant features at 94, 131 and 335 A after 
the main impulsive peak of the flare. The bright source located between the ribbons, 
visible at the impulsive phase m_all AIA wavelengths and also seen in HXR (Section 
0]), was analysed in detail bv lSimoes. Graham, and Fletcher ( 2015 ). They characterised 
this source as the main energy release site of this event. We summarise their findings as 
follows: 1) the source is located in the corona (possibly low down), it has a filamentary 
shape along the loops seen in EUV/UV images, the lack of hot loops connecting the 
region after the impulsive phase, and the weak and featureless photospheric magnetic 
field at the same location (see Figure 0]). 2) it has intense and impulsive EUV and HXR 
emission. 3) consistently, the source is found to be dense (log n- 11.50 ± 0.82) and hot 
(T « 12 ~ 16 MK). 4) strong red-shifts observed in many EUV emission lines observed 
by Hinodel EIS, including Fe xn and Fe xxiv, indicate plasma downflows of 40-250 km 
s" 1 , which, along with plasma outflows observed in AIA images, are interpreted as 
plasma outflows along the magnetic loops. 

We investigate the emission from the ribbons and the coronal source (for context) 
separately by marking regions of interest (ROI) associated with each of the three main 
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Figure 2. SDO/AIA images at 94, 131, 171 and 1600 A filters of SOL2013-11-09, at four times during the 
impulsive phase. 


sources as shown in Figure [5J East (orange) and West (green) ribbons, coronal source 
(blue). We obtained lightcurves for all AIA passbands by summing up the pixels for 
each ROI, which are shown in Figure [6j We also show RHESSI HXR 15-25 keV 
counts for comparison. The UV (1600 and 1700 A), 304 A, and the “warm” EUV 
(171, 193, 211 A) channels have similar enhancements during the impulsive phase, 
showing impulsive emission from both ribbons and the coronal source. The bursts are 
coincident with the HXR main peak (06:25:46 UT), within the instrumental cadence. A 
secondary, less impulsive HXR peak at 06:27:14 UT is also well-associated with most 
AIA channels. The peaks are simultaneous at the three ROIs in all nine AIA filters, 
within the instrumental cadence. After about 06:30 UT the coronal source fades out at 
all wavelengths. At 94 and 131 A the coronal source ROI is dominated by the emission 
from the coronal loops. 

3.2. Emission Measure Distributions 


To investigate the evolution of the plasma at the flaring ribbons^we apply a method of 
regularised inversion to the AIA data developed bv [Hannah and Kontarl ( 2012 ) to obtain 
the differential emission measure (DEM), £(T) = n e 2 dh/dT [cm -5 K _1 ], defining the 
quantity of emitting material as a function of temperature, T , along the given line-of- 
sight, h, with an average d ensit y n e . We employed up-to-date AIA temperature response 
functions ( Boerner etal\ 120121. 2014 ) which have e mpirical corrections for the missin g 
emission lines from the CHIANTI v7.1.3 database dDere et al, . 1 997l : lLandi et <3/i2013l) . 
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Figure 3. SDO/AIA images at 193, 211, 335 and 304 A filters of SOL2013-11-09, at four times during the 
impulsive phase. 


time-dependent response corrections for each channel due to degradation of the filters, 
and normalisation to ensure agreement with SDO/EVE spectroscopic full-disk obser¬ 
vations. The regularisation method also provides the DEM uncertainty and effective 
temperature resolution. To obtain the DEM it must be assumed that the emitting plasma 
is optically thin, in local thermodynamic equilibrium, and in ionisation equilibrium. In 
the case of a flare ribbon, where the emitting plasma is both hot and dense, this may 
be hold, although departures fro m equilib ri um and optical depth effects should never 
be discounted completely (see lGraham et al\\20 1 3l. fo r further discussion). C ontinuum 
emission may also contribute to the AIA passbands dO’Dwver et al. lEoioh . however 
this contribution is likely to b e small for C class flar es and can probably be neg lected 
( Milligan and McElrovll2013 ). We point the reader to lHannah and Kontarl ( 2012 ) for an 
extended explanation of the method and comparison with other available methods. 

We applied the method to the average emission inside each ROI indicated in Figure[5] 
throughout the impulsive phase. Integrating the DEM over fixed temperature intervals 
gives the emission measure distribution, EMD = n 2 h, in the more practical units of 
cm -5 ( e.g . iGraham et a/.ll2013l : lHannah and Kontadl2013l) . The resulting EMD for the 
East and West ribbons and coronal source for selected times are shown in Figure[7]along 
with the EM -loci curves for each filter, these representing the EMD for an isothermal 
plasma, i.e the maximum theoretical EM. 

The pre-flare EMD is shown by the dashed line with a peak around log T = 6.3 ~ 6.4 
and falling off towards higher temperatures. A second peak is visible around log T = 
6.9 ~ 7.1, most notable in the East ribbon, however the 94 and 131 A AIA channels 
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Figure 4. SDO/HMI line-of-sight (LOS) magnetogram, overlaid by SDO/A IA 1 6 00 A contours a t 500 DN 
s -1 , showing the position of the ribbons and the coronal source, analysed bv ISimoes. Graham, and Fletched 
(2015). 
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Figure 5. a) Regions of interest (ROI) defined for the AIA images: East (orange) and West (green) ribbons, 
and coronal source (blue). 
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Figure 6. SDO/AIA lightcurves for EUV/UV filters, (a) 94 A, (b) 131 A, (c) 171 A and (d) 1600 A, for each 
region of interest (ROI), as defined in Figure \5\ East (orange) and West (green) ribbons, and coronal source 
(blue). RHESSI 15-25 keV counts (grey) are shown as a reference of the impulsive non-thermal emission, 
with the time of its maximum shown by the vertical dotted line. 


both contain contributions from lines at approximately 1 MK and 10 MK, therefore it 
is likely that this apparent hot emission is in fact much cooler. The EMD for the three 
sources share a similar development; a pre-flare EMD which evolves with an overall 
increase of the emission measure at all temperatures in the AIA passband response 
range (roughly 5.7 < log T < 7.3), a stronger increase in high temperature emission at 
log T - 7.0, and the formation of a low temperature “shoulder” just below log T = 6.0. 

We then integrated the EMD at temperature ranges 0.5-1.5 MK, 1.8-3.2 MK and 
7.9-12.6 MK, which show the most prominent peaks, to obtain the column emission 
measure, EM C = n 2 h [cm -5 ], at these temperature ranges. We did not attempt this at 
the temperatures around log T = 6.6, as this range is not very well constrained by the 
AIA filters. In order to compare these results with (volume), EM = n 2 V [cm -3 ], values 
from GOES and RHESSI, we converted the EM C to EM by multiplying the values by 
the projected area, A, of the sources, i.e. making V - hA. The projected area A of 
the emitting plasma was estimated by finding the number of pixels inside each ROI 
with a value above a determined threshold. Here we used 94 A images and chose a 
threshold value of 35 DN s -1 pixel -1 that captures the weaker but relevant emission at 
the early phase, and also the stronger emission at the impulsive phase. We verified if 
the emitting pixels at the other AIA filters sensitive to cooler plasma temperatures co¬ 
aligned well with the emitting area at 94 A, generally finding a good agreement within 
a couple of pixels. The time evolution of the EM (now in cm -3 ) at the three temperature 
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Figure 7. Emission measure distribution EMD (cm -5 ) for the East (orange) and West (green) ribbons and 
coronal source (blue) at selected times from the early phase into the impulsive phase. Pre-flare (06:20:19 UT) 
EMD are shown by the grey dashed line for comparison. The EM -loci curves for each filter are indicated by 
different colours and identified in the top-left panel. 


ranges for the three sources is shown in Figure [8] along with the values from GOES 
and RHESSI obtained for full-Sun observations. We subtracted the EM pre-flare values 
to consider the flare excess only, taking the pre-flare values as the mean in the interval 
6:00-6:20 UT. Looking at Figures [S^-c in the early phase (6:22:00-6:25:20 UT), the 
EM for the two ribbons in all three temperature bands starts to rise very slowly from 
about 6:22 UT (clearly visible on a logarithmic scale, not shown here), with a faster 
enhancement after 6:24:40 UT. During the main impulsive phase (6:25:20-6:26:40 UT) 
the EM in the three temperature ranges for all the three sources rises impulsively, 
peaking around 6:25:40-6:26:20 UT. The GOES and RHESSI emission measures also 
show a steeper enhancement. This phase coincides with the main HXR peak. The rapid 
increase, followed by a decrease of the EM implies fast heating and cooling (or removal) 
of the plasma, and it is more pronounced at the East ribbon. After 6:26:40 UT, at 7.9- 
12.6 MK, the EM keeps rising, which is consistent with loop filling by chromospheric 
evaporation. AIA images at 94, 131 and 335 A show the enhancement of bright loops 
in this phase, and the increase is also supported by the GOES and RHESSI EM values. 

At lower temperature ranges, the EM peak is more pronounced, showing a fast in¬ 
crease and decrease in the amount of material below « 3 MK. This behavior is similar 
to emission lines observed by SDO/EVE formed at chromospheric and transition region 
temperatures, 4.2 < log T < 5.7K, which show impulsive emission in time with HXR. 
In Figure[8]i we show the lightcurves of C hi 977 A, O iv 554 A and O v 629 A, noting 
that other lines formed at similar temperatures display similar behaviour, namely Ne vn 
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Figure 8. Time evolution of the EM (converted to cm -3 for comparison with GOES and RHESSI results) at 
three temperature ranges, (a) 0.5-1.4, (b) 1.8-3.2, and (c) 7.9-12.6 MK, for East ribbon, West ribbon, coronal 
source from AIA from pre-flare into the impulsive phase, (d) SDO/EVE lightcurves of selected emission lines 
formed at transition region temperature (log T « 5). The vertical dotted line indicates the time of the HXR 
peak. 


465 A, O hi 526 A, He i 584 A, O hi 599 A, O iv790 A, Q vi 1032 A, w ith s ome being 
noisier and less impulsive. Simulations by [Fisher. Canfield, and McClvmontl ( 1985 ) 
show that the entire flare chromosphere cools rapidly on a radiative timescale after the 
heating is turned off. However, in this flare, after the main HXR peak there is still HXR 
and EUV emission present, indicating that the energy deposition is still occurring. 


4. RHESSI Hard X-rays Observations 


We now employ RHESSI HXR observations to characterise the hot plasma and the non- 
thermal electrons. First, we evaluate the HXR spectra during the impulsive phase, with 
a 4 second time resolution. RHESSI spectra were fitted with an isothermal plus a sin¬ 
gle po wer-law thick-target model, using Object Spectral Executive (OSPEX) software 
( Schwartz et al. 20021) . No attenuator was active during the course of the event (state 
AO). The non-thermal power-law is assumed to originate from a cold, collisional thick- 
target ferownlll971 ). The fitting parameters for the non-thermal electrons derived from 
RHESSI are shown in Figure [9] We will use these parameters to derive the collisional 
beam heating in Section [5^21 The thermal plasma parameters (EM and T) are shown in 
Figure \T\ along with the same parameters estimated from GOES data. 
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Figure 9. Spatially unresolved RHESSI spectral analysis results for the non-thermal electrons: F tot , the total 
electron rate above the low energy cutoff E c and spectral index S. 


We applied the imaging spectroscopy technique for the three main HXR sources 
observed associated with the EUV ribbons and coronal source, indicated in Figure [TO] 
Due to RHESSI dynamic range, it was not possible to do imaging spectroscopy for the 
West ribbon. We constructed RHESSI Clean images for 14 energy bands, logarithmi¬ 
cally spaced between 3 and 40 keV, using detectors 3 to 8, integrated during the main 
impulsive phase 06:23:34 - 06:26:42 UT. The spectra were fitted with an isothermal 
plus thick-target model, shown in Figure \Tl\ and the resulting parameters are found 
in Table [T] For comparison, we also fitted the spatially unresolved full Sun spectrum 
integrated for the same time interval. The thermal and non-thermal properties of the 
three HXR sources are very similar, confirming that the ribbon plasma is above 10 
MK during the impulsive phase, as indicated by our AIA analysis in Section [3711 Also, 
the non -thermal properties are in agreement with the findings of Simoes. Graham, and 
Fletcher ( 2015h who showed that the coronal source is 9-18 arcsecs long, with a plasma 
density of log n- 11.5 and thus being collisionally thick for electrons with energies of 
up to 45-65 keV. 


5. Discussion 

5.1. Impulsive Heating of Ribbons 

We now investigate the evolution of the properties of the thermal plasma at the East 
ribbon and coronal source compared to the full Sun spectral results. The presence of 
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Figure 10. (a) RHESSI HXR 23-27 keV contours overlaid on an SDO/AIA 171 A image (inverted colours) 
with (red) dashed boxes indicating the regions for imaging spectroscopy, (b) RHESSI HXR 4-9 keV contours 
obtained with Clean (red) and MEM NJIT (dark blue) overlaid on an SDO/AIA 131 A image (inverted 
colours). 



photon energy [keV] photon energy [keV] 

Figure 11. RHESSI photon spectra for the three main HXR sources indicated in Figure Q0| and spatially 
integrated, for the interval 06:23:34 - 06:26:42 UT. 


Table 1. RHESSI imaging spectroscopy results. 



EM xl0 47 cm -3 

T[ MK] 

F 35 IO 35 s -1 

5 

East ribbon 

1.0 ±0.5 

10 + 1 

1.0 + 0.3 

5.7+ 0.3 

Coronal 

0.6 ± 0.3 

12 + 1 

0.7 + 0.2 

5.2 + 0.3 

West ribbon 

1.0 ±0.7 

10 + 1 

0.5 + 0.2 

4.9 + 0.3 

Full Sun 

1.33 ±0.02 

11.81 ±0.05 

2.84 + 0.04 

5.61 + 0.02 
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low-energy HXR from thermal bremsstrahlung from the ribbon sources is shown by 
imaging spectroscopy results, and can also be seen directly in the reconstructed images. 
HXR 4-9 keV emission associated with both ribbons (and coronal source) is shown 
in Figure [TOb. We obtained images using two different imaging algorithms, Clean and 
the Maximum Entropy Method of the New Jersey Institute of Technology (MEM-NJIT: 
ISchmahl et al, ■ 2007 ), to confirm that the sources were not spurious artifacts of the 
algorithms. MEM-NJIT has a tendency to over resolve the sources, but in this case it 
works to confirm the overall spatial distribution of the low-energy HXR emission over 
part of the ribbons. We then constructed Clean maps (with Clean_beam_widtl{3 =1.5; 
detectors 3-8) for nine time bins in the interval 06:22:22 and 06:27:51 UT. The spectra 
were fitted with an isothermal plus thick-target component. The latter was included 
to achieve a better fitting for the temperature, as fluxes around 8-10 keV and above 
could not be fitted with a thermal component only, requiring a non-thermal tail. In 
Figure [12] we show the EM and T for the East ribbon (orange), coronal source (blue) 
and full Sun (grey). Although the uncertainties in the imaging spectroscopy method 
are larger due to the smaller time intervals considered, the temperature peak around 
06:25:40 UT is evident for both East ribbon and coronal source, in good agreement with 
the full Sun values from both RHESSI and GOES (see FigureQil). The evolution of the 
plasma temperature in the impulsive phase of this event is in fact quite unusual. The 
plasma temperature from GOES and RHESSI (Figure |T4]) show that the temperature 
peak (06:25:38 UT) occurs about 8 to 10 second^ before the HXR peak (06:25:46 UT). 
Even considering the different instrumental responses from GOES and RHESSI, both 
sample the hottest plasma present in flares, weighted by the EM jRvan et all 20140 , 
indicating that this temperature peak reveals an impulsive heating followed by a fast 
cooling process. The temperature decrease from 13 MK to 11MK in about two minutes 
is consistent with conduction l osses, w ith a time sc ale of about 14 to 38 seconds (for 
the coronal source. iSimoes . Graham , and Fletchen 120151) . The e-folding cooling time 
for the (full Sun) temperatures in Figure [14] is t - 35 - 55 seconds. 

Since the spatial resolution of RHESSI is lower than that of AIA, we confirm the 
location of the hi ghest EM by applying the AIA-DEM method for each pixel in AIA 
images, following lHannah and Kontan d2013l) . In Figure [l3l we show EM maps from 
the AIA-DEM method, where the largest concentrations (EM C > 10 29,5 cm -5 ) of hot 
plasma (T > 8 MK) are near the flare ribbons and the coronal source. 

EM maps (Figure[l3]) obtained with the AIA-DEM method also suggest the presence 
of short-lived, high-EM ribbon (EM C « 10 30 cm -5 ) at 11-14 MK near the HXR peak 
time at 06:25:46 UT. 

The gradual rise of the EM values from the AIA-DEM (Figure [8]) starting about 
06:24 UT shows the filling of the loops with hot material but the sudden enhancement 
between 06:25:46 and 06:25:55 UT comes from this impulsively heated material, with 
the two components then contributing to the overall EM. The fast cooling of the 13 
MK source can only be attributed to conduction to the lower atmosphere, as radiative 
cooling is much slower than the observed e-folding time, of about 35-55 seconds. It 


1 Arbitrary factor applied to the beam convolved with Clean sour ces. Values between 1.5- 2.0 tend to give 
better results when compared to different imaging algorithms (see ISimoes and Kontadl2013[ and references 
therein.) 

2 Given the time resolution of GOES (« 2 s) and RHESSI (« 4 s) 
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Figure 12. EM and T from RHESSI imaging spectroscopy for the coronal source (blue crosses) and East 
ribbon source (orange). The EM and T from the unresolved spectral fitting are also shown (grey). 


is not clear what is the source of energy to produce the temperature peak, however we 
here speculate a scenario where a rise in the plasma temperature leads the acceleration 
of particles. 

5.2. Energy Budget at the 10 MK Ribbon 


We now examine the energy budget of the 10 MK plasma in the East ribbon. For that, we 
consider a slab of plasma at T = 10 MK in the lower atmosphere (i.e. at transition region 
heights), with a thickness L, being collisionally heated by the non-thermal electrons, 
balanced by radiative and conductive losses. 

Following Fletcher et al . ( 20131) we calculate the energy budget of the 10 MK plasma, | 
considering radiative and conductive losses and the collisional energy loss of the non- 
thermal electrons as the energy input mechanism. The total collisional power input 
Peon [erg s -1 ] to the plasm a can be inferred from the collision al thick-target model by 
integrating the energy loss dEmslielll978l : lFletcher et al, Ibonh 


Peon — FtotEc 


6-1 


6-2 




1 - 6/2 

x c 



( 1 ) 


where B is the incomplete beta function and v c = 3KN/E 2 , where K = 2ne 4 , and also 


F tot = f FoE^dEo (2) 

Je c 

being E tot the total number of electrons per second above E c injected into the source. 

The values for E tot , 6 and E c are taken from RHESSI spectral analysis (Figure©. In 
order to estimate the time evolution of these parameters for the East ribbon source, we 
used the results from the full Sun spectroscopy and assumed that about a third of the 
total number of electrons are associated with each of the three HXR sources, based on 
the HXR imaging spectroscopy results in Table [T] 

The column depth, N, can be estimated from the column, EM C , inferred with the 
AIA-DEM analysis, as N = n e L - (EM L) 1/2 , for a uniform source of thickness L. 
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Figure 13. EM maps obtained from the AIA-DEM method dHannah and Kontar]l2012l) . for temperature 
ranges 2-3, 6-8, 8-11 and 11-14 MK, at four times during the impulsive phase. The arrows (from left to 
right) in panel k indicate the positions of the East ribbon, coronal source and West ribbon. 


The hot plasma will lose energy by conduction, along the field direction z, to cooler 
layers of the atmosphere at a rate L con d [erg^s -1 ]. We estimate the conductive losses 
following [Battaglia. Fletcher, and Bend ( 2009 ). considering the case of flux-limited 
conduction: 


r c0 nd = q(x)k 0 T 5,2 ^A 
dz 


(3) 


where the factor g(x), a function of x = log(/ m f p /L), reduces the classical Spitzer 19651) 


conduction coefficient kq = 


10 6 erg cm 1 


,-i 


K 7/2 to the flux-limited conduction 


dCampbell 19841) . Battaglia. Fletcher, and Benzl ( 2009 ) fitted the values of g published 
by Campbell! ( 19841) as g(x) = 1.01e _005( * +6 - 63)2 . This condition applies because the 
chromospheric scale length is small and the electron collisional mean free path, / m f p = 
5.21 x 10 3 T 2 /n e , is significant compared to the temperature scale length. Here, we 
approximate dT I dz by T/L and set T = 10 MK and L - 2000 km. 

The optically thin hot plasma ra diates energy at a rate L rad that can be estimated by 
( Rosner. Tucker, and Vaianall 1978 ): 


T rad = 10“ 17 - 73 EM T~ 2 '\ (10 6 - 3 < T < 10 7 K). 


(4) 


Using the EM values derived from the AIA-DEM analysis (Figure [8]:), and the area 
A of the 10 MK source estimated as defined in Section [3711 we found that the collisional 
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Start Time (09-Nov-13 06:24:00) 

Figure 14. Plasma temperature inferred from RHESSI (blue) and GOES (orange) data, peaking about 10 
seconds before the RHESSI HXR counts at 3-6 keV (green) and 25-45 keV (violet). The values of the 
e-folding cooling time, r (in seconds), are indicated in the figure. 


heating is not sufficient to balance the estimated conductive losses as shown in Figure 

EU _ 

Fletcher et al. ( 2013 ) found a substantial amount of plasma at 10 MK in the flare 
ribbons during the early impulsive phase of the flare SOL2010-08-07T18:24, with an 
average column EM of a few times 10 28 cm -5 . They found that the energy carried by an 
electron beam is not sufficient to heat the ribbon plasma when radiative and conductive 
losses are taken into account, unless a low-energy cutoff of E c « 5 keV is considered. 
The total electron energy inferred from the HXR spect rum is mostly defined by the 
low-energy cutoff E c (e.g. Saint-Hilaire and Benz 120051) . which is sometimes chosen 
arbitrarily to balance the energy deposited by th e non- thermal electrons and the value 
of the maximum thermal energy (e.g. iMrozek and Tomczakll2004l) . The E c values we 
found from the HXR spectral fitting are already quite low, approximately 8 keV during 
the impulsive phase (see the bottom panel in Figure O, only a few times the average 
thermal energy kT of electrons in the 10 MK plasma. It is of course possible to reduce 
the value of the low-energy cutoff to make the energy losses and gains balance, and a 
lower value of the low-energy cutoff is permitted by the HXR spectrum. A lower E c 
means an increased total ‘non-thermal’ power, and an increased electron number flux. 
For this event, a value of E c « 4 keV is necessary to equate energy gains and losses. 
On the other hand, using the E c values obtained from the HXR spectral analysis, a slab 
thickness L of ^30 Mm is required to balance the energy budget, which would comprise 
a large portion of the loops, a picture not supported by the AIA images, which show 
very compact footpoint sources. 
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Figure 15. Power budget of the East ribbon source, considering beam collisional heating (blue), conductive 
(magenta) and radiative (gold) cooling, in a Lj - 2000 km thick slab of plasma at a temperature of 10 MK. 


As shown bylSimoes, Graham, and Fletcher (2015), the coronal source is both dense 
(n « 10 11 cm -3 ) and hot (T « 13 MK), a downward thermal conduction from the coronal 
sour ce to the ribbons may contribute to the observe d heating at thes e regions (e. 2 . Hori 
et fl/. ll998l : l0iu et <3/.ll2013l : [Battaglia, Fletcher, and Benzll2009h . 

If the low energy cutoff is indeed as low as 4 keV to balance energy gains and losses 
by the chromospheric plasma then E c approaches the mean thermal energy of the 10 
MK plasma, and it becomes difficult to separate the ‘beam’ from the thermal core. We 
may find it profitable in the future to think about the electrons in the chromosphere 
forming a single distribution, such as a /^-distribution, which describes a single pop¬ 
ulation transitioning smoothly from a power-law-like tail to a thermal-like core. This 
distribution is being discussed in connection with dense coronal HXR sources, where 
populations of heated ele c trons and ac c elerate d electrons are present in the same volume 
( Oka et JZll2013l. [2015 ). Bian et al. ( 2014 ) have demonstrated how a ^-distribution 
arises naturally in a volume where diffusive acceleration and collisional losses operate 
simultaneously and co-spatially. If this volume is coronal and a beam is formed by a 
high energy tail that ‘leaks out’ then that beam looks very much like a standard power- 
law beam, with a low energy cutoff determined by the particle velocity at which the 
escape time is less than the acceleration or diffusion timescales; then, we are then still 
in the standard collisional thick-target scenario but with a slightly differently shaped 
coronal beam. However, if we move the site of both heating and acceleration to the 
chromosphere, which is a very different interpretation, then the properties of the Re¬ 
distribution extending across all of the electron energy space can be investigated to 
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evaluate the total required power and number flux needed to account for all radiation 
signatures. This is an interesting possibility; it does however require that some other 
agent, such as wave turbulence, is present in the chromosphere to locally heat and 
accelerate electrons. A full investigation is beyond the scope of this paper. 


6. Summary 


We have presented an analysis of the plasma in the flare ribbons of the event SOL2013- 
11-09, a C2.6 class non-eruptive, two-ribbon flare using SDO/AIA EUV and RHESSI 
HXR observations. The ribbons have impulsive EUV/UV emission seen in all SDO/AIA§ 
filters, well associated with non-thermal HXR emissio n obs erved by_RHESSI. Using 
the method of regularised inversion of SDO/AIA data faannah and Kontanl2012l) we 
obtained the differential emission measure (DEM) of the two flare ribbons, and investi¬ 
gated the time evolution of the emission measure (EM) in three temperature ranges (0.5 
- 1.4, 1.8 - 3.2 and 7.9 - 12.6 MK). From these, we have shown that the plasma heats 
rapidly to 12-13 MK during the impulsive phase of the event, marked by the HXR peak. 
The EM temporal evolution shows a peak near the time of maximum HXR emission, in¬ 
dicating fast heating and cooling with the hottest plasma (T =7.9- 12.6 MK) reaching 
EM values of 1 - 3 xlO 47 cm -3 , these values agreeing with those obtained from GOES 
and RHESSI. The rapid evolution of the ribbon plasma temperature and its high peak 
temperature are confirmed by RHESSI imaging spectroscopic analysis, and also agree 
with the temperature derived from GOES (although without spatial resolution). Also 
we note that the evolution of the ribbon plasma characteristics are very similar to the 
those found in the int ense and compact coronal source, as studied in detail bv Simoes. 
Graham, and Fletcher d2015l) . Performing RHESSI HXR imaging spectroscopy, we 
obtained the parameters to describe the distribution of non-thermal electrons at each 
source (both ribbons and coronal source). With the information about the plasma and 
non-thermal electrons at the East ribb on, we deduced the energy balance o f the plasma, 
considering collisional beam heating dEmsliel 1 197 8l : iFletcher et r//.ll2013l) against con¬ 
ductive and radiative losses. We found that beam heating alone is not sufficient to heat 
and maintain the ribbon plasma at T = 10 MK, even with a low energy cutoff of E c « 8 
keV, and speculate if the de nse and hot coronal source can pr ovide a heat source for the 
ribbons by conduction ( e.g . iHori et al 1998: Pin et al H 2013 I) . 
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